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We outline the expected constraints on non-Gaussianity from the cosmic microwave background 
(CMB) with current and future experiments, focusing on both the third (/nl) and fourth-order 
(gNL and tnl) amplitudes of the local configuration or non-Gaussianity. The experimental focus is 
the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted 
maps. In adition to a measurement of tnl and gtih with WMAP 5-year data, our study provides 
the first forecasts for future constraints on (7nl. We describe how these statistics can be corrected 
for the mask and cut-sky through a window function, bypassing the need to compute linear terms 
that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness 
estimator. We discus the ratio Anl = 'nML/(6/NL/5)^ as an additional test of single-field inflationary 
models and discuss the physical significance of each statistic. Using these estimators with WMAP 
5- Year V-|-W-band data out to Zmax = 600 we constrain the cubic order non-Gaussianity parameters 
TNL, and gNL and find —7.4 < ^nl/IO^ < 8.2 and —0.6 < tnl/IO* < 3.3 improving the previous 
COBE-based limit on tnl < 10* nearly four orders of magnitude with WMAP. 



I. INTRODUCTION 

We have now entered an exciting time in cosmological 
studies where we are now beginning to constrain simple 
slow-roll inflationary models with high precision obser- 
vations of the cosmic microwave background (CMB) and 
large-scale structure. In addition to constraining infla- 
tionary model parameter space with traditional parame- 
ters such as the spectral index and the tensor-to-scalar 
ratio r, we may soon be able to use parameters associ- 
ated with primordial non-Gaussianity to improve model 
selection. 

In the simplest realistic inflationary models, the field(s) 
responsible for inflation have minimal interactions. Such 
an interaction-less situation should have led to Gaussian 
primordial curvature perturbations, assuming that pertn- 
bations in the inflaton field generates the curvature per- 
turbation. In this case, the two point correlation function 
contains all the informations on these perturbations. If 
the early inflation field(s) have non-trivial interactions, 
higher-order correlation functions of the curvature per- 
turbations will contain connected pieces encoding in- 
formation about the primordial inflationary interactions. 
This is analogous to the situation encountered in particle 
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physics where correlation functions can be separated into 
unconnected and connected Feynman diagrams, the later 
containing information about the underlying interactions 
(see Fig. [T] for an example involving the four-point func- 
tion). A detection of non-Gaussianity therefore gives an 
important window into the nature of the inflation field(s) 
and their interactions. 

To parameterize the non-Gaussianity of a nearly Gaus- 
sian fleld, such as the primordial curvature perturbations 
C(x), we can expand them perturbatively to second 
order as: 

C(x) = Cs(x) + [C,'(x) - (C,^(x))] + ^5nlC,'(x), 

(1) 

where Cg(x) is the purely Gaussian part with /nl and 
.gNL parametrizing the first and second order deviations 
from Gaussianity. This parameterization of the curva- 
ture perturbations is known as the local model as this 
definition is local in space. 

Much effort has already gone into measuring non- 
Gaussianity at first-order in curvature perturbations us- 
ing the bispectrum of the CMB anisotropies or large-scale 
structure galaxy distribution parametrerized by /nl (see 
Eq. [T]). These studies have found /nl to be consistent 
with zero [IHlj- However, there is hope that a significant 
detection may be possible by future surveys that will lead 
to improved errors 5J. 

In the trispectrum, two parameters of second-order 
non-Gaussianity at fourth-order in curvature perturba- 
tions, Tnl and ^nl, can be measured. In this paper 
we also introduce a third parameter, Anl is an addi- 
tional parameter that compares tnl of the trispectrum 
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ical significance of eacli statistic. In Section III 



describe the skewness and kurtosis power spectra and 
explain how they may be used to extract information 
about primordial non-Gaussianity from the CMB. In Sec- 



tion IV, we describe the signal-to-noise of each estimator, 
how to add the experimental beam and noise to these cal- 
culations and discuss why these power spectra have the 
advantage for dealing with a cut sky. In Section [V] we 
calculate the fisher bounds for upcoming experiments for 
each statistic. In Section IVTl we discuss the technical de- 
tails for measuring non-Gaussianity in the trispectrum 
and in Section [VIII we conclude with a discussion. 



FIG. 1: Four point correlation function for the (f>^ theory. 
The correlation functions breaks up into interaction-less un- 
connected diagrams and connected diagrams containing in- 
formation about the interactions. 



to (6/nl/5)^ from the bispectrum as a ratio: 

(6/nl/5)^ 

This ratio can be quite different for many inflationary 
models [6l[7] and, as will be shown below, Anl 7^ 1 rules 
out single-field inflationary models altogether, including 
the standard curvaton scenario (which neglects pertur- 
bations from the inflaton field). 

In this paper we discuss the skewness and kurto- 
sis power spectra method for probing primordial non- 
Gaussianity and give constraints for the first (/nl) and 
second-order (^nl and tnl) amplitudes of the local model 
in addition to their ratio Anl- Using the bispectrum of 
CMB anisotropics as seen by WMAP 5-year data, Smidt 
et al. (2009) found -36.4 < /nl < 58.4 at 95% confi- 
dence This is to be compared with the most recent 
WMAP 7 measurement of -10 < /nl < 74 [3], where 
part of the discrepancy is due to a difference in optimiza- 
tion [8]. As outlined in Section VI, using the trispectrum 
of the same data we find that —0.6 < tnl/IO"* < 3.3 
and -7.4 < .gNL/10^ < 8.2 at 95% confidence level show- 
ing second order non-Gaussianity is consistent with zero 
in WMAP. This paper serves as a guide to the analysis 
process behind our derived limits on tnl, <7nl and ^nl- 

Furthermore, in this paper we analyze what to real- 
istically expect when measuring non-Gaussianity from 
CMB temperature data. We believe establishing what 
constraints can be placed upon /nl, tnl, 3nl and Anl 
by future experiments is important in determining what 
models may and may not be tested by future data. We 
also highlight several advantages of our work, including 
ways to correct the cut-sky and mask through a window 
function without using linear terms which are computa- 
tionally prohibative [9, 10 . 

This paper is organized as follows: In Section |ll] we 
review how non-Gaussianity may be used to distinguish 
between common inflationary models and stress the phys- 



II. NON-GAUSSIANITY FROM COMMON 
INFLATIONARY MODELS 

Non-Gaussinity is a powerful tool that may be used 
to distinguish between inflationary models. The sim- 
plest models do not produce a detectable amount of non- 
Gaussianity. Maldacena jTTj has shown that a single- 
field, experiencing slow roll with canonical kinetic energy 
and an initial Bunch-Davies vacuum state produces 

hL^^{ns + f{k)nt). (3) 

Here rig and nt are the scalar and tensor spectral indices 
respectively. The function f{k) has a range < f{k) < | 
based on the triangle shapes (see below) of the fci such 
that / = in the squeezed limit and / = | for an equilat- 
eral triangle. For this reason, /nl < 1 will remain unde- 
tectable in the simple slow roll scenario with CMB data 
alone. If any of the above assumptions are violated, very 
specific types of non-Gaussianity are produced [5l [T2l[T3] . 
In the bispectrum B(^(ki, fc2, ^3) defined by 

(CkiCk.Ck3) = (27r)35(ki + k2 + k3)i3c(fci, A:2, fca), (4) 

where C is the primordial curvature perturbation, non- 
Gaussianities show up as triangles in Fourier space. Dif- 
ferent triangle shapes are be produced by different un- 
derlying physics, for example: 

• squeezed triangle (fci ^ fc2 3> /ca) This is the domi- 
nating shape from multi-field, curvaton, inhomoge- 
neous reheating and Ekpyrotic models. 

• equilateral triangle (fci = fc2 = k^) This 
shape is produced by non-canonical kinetic energy 
with higher derivative interactions and non-trivial 
speeds of sound. 

• folded triangle {ki = 2^2 = 2^3) These triangles are 
produced by non-adiabatic-vacuum models. 

Additionally, linear combinations of the above shapes or 
intermediate cases such as elongated triangles (fci = k2 + 
k^) and isosceles triangles (fci > ^2 = ^3) are possible [SI 
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FIG. 2: Plot of the shape functions 5''°'^^'(1, fe, ^3) and 
5'°''""(l,fc2,fc3) normahzed such that ^(l, 1, 1) = 1. In these 
plots only values satisfying the triangle inequality ^2 + fcs > 
fci = 1 as well as the requirement k2 < fca to prevent show- 
ing equivalent configurations are non-zero. The plot on top 
verified 5'°='^' is maximized when fci ~ fca ^ k2 whereas the 
bottom plot verifies 5"='?"''' is maximised when fci ~ fc2 ~ ^3. 



[T2I [T3] . The most recent WMAP 7 constraints on the 
amount of non-Gaussinaity from each shape is —10 < 



at 95% confidence ^ 



214 < /nl 
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< 266 and -410 < f\ 



orthog 
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common configurations as: 
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Km' 



(6) 
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5™(fci,fe,fc3) (X ^(K,2-Ks)+4J^,{8) 



Ki 



11 



where 



Kr, 



K, 



pq 



^{k,)P with K = Ki, 

i 



(9) 
(10) 

(11) 



hp = Kp-2{h)P with h^ka, (12) 



with Ar 



1 + Spq and Apqr 



Apqi^Aqrp 



(no 



summation). Plots for the local and equilateral shapes 
are given in Figure [2] 

In addition to /nl being generated by different shapes, 
it also may vary with scale. Recently, a new parameter 
has been introduced to measure this scale dependance 
defined as: 



rfln|/NL(fc)| 
din fc 



(13) 



This scale dependance has the ability to test the ansatz[T] 
to test whether the local model should allow for /nl 
to vary with scale [15]. Using the results of Smidt et 
al.(2009) (Fig. 16 of Ref [5) and assuming 
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(14) 



we can constrain nf^^{l) to roughly —2.5 < nf^^{l) < 2.3 
at 95% confidence. We therefore find /nl is consistent 
with having no scale dependance. 

In this paper we focus on the local model that probes 
non-Gaussianty of a squeezed shape. As mentioned 
above, simple inflationary models can not produce a de- 
tectable amount of non-Gaussinity for local models. We 
now review the prediction for local non-Gaussianity for 
the most common models. 



A convenient way to distinguish between shapes is to 
introduce the shape function defined as 

S{ki,k2,k3) = -^(fcifc2fc3)^S<j(fci,fc2,fc3), (5) 

where N is a normalization factor often taken to be 
1//nl- Using a notation introduced by Fergusson and 
Shellard |14| , we can give the shape function for the more 



A. Review Of The SN formalism. 

The curvature perturbation can be conveniently de- 
scribed using the SN formalism [T6Vl20j . During infla- 
tion, spacetime expands by a certain number of e-folds 
N. By Heisenberg's uncertainty principle, expansion for 
each point in space ends at slightly different times pro- 
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ducing a spatially dependent total e-fold: 



N{x) = J H{x,t)dt, 



(15) 



where H{x, t) is the Hubble parameter allowing us to de- 
fine N{x) — N + SN{x). The fluctuations in e-fold about 
the mean value N, which correspond to perturbations in 
local expansion, are the curvature perturbations C — ^N. 

In addition to a spatial parameterization, we may pa- 
rameterize the number of e-folds by the underlying fields 
C = N{(j)^) — N where (j)^ represents the initial val- 
ues for the scalar fields . If we write out the fields as 
(j)^ = (j)^ + Scj)"^ we can expand the curvature perturba- 
tions as 



(16) 



The Nx means the derivative of N with respect to the 
fields X. For example, NaiA2 = q^a^ q^a2 ■ In this equa- 
tion there is an implicit sum over the Ai. Einstein sum- 
mation is implicit in all equations relating to the 6N 
formalism. 

Using this formalism we may compute to first order 
from C — NaSip^: 

(CkCk') = NaNbC^^ {k){2nf 6' {k + k'), (17) 

where C^^{k) in the slow roll limit becomes to leading 
order S^^P{k). 

Likewise, we can calculate the bispectrum and trispec- 
trum in this formalism; 

Bc{ki,k2,k3) = NANBcND[C^''{ki)C''''{k2) (18) 

+ C^^(fci)C^^(fc2) + C^^(fci)C^^(fc2)], 

T^{ki,k2,k3,k4) = Na,a,Nb,b,NcNd (19) 
x[C^'^Hki3)C^'^{k2)C^'°ik2) + (11 perms)] 

+Na,a2A,NbNcNd 
X [C^'''{ki3)C^'''{k2)C^'''{k2) + (3 perms)], 

where fc^ = jk^ -I- kj|. In the slow roll limit to leading 
order these expressions may be rewritten as: 



values for each statistic: 

5 NaNbN"^^ 



/nl — 



5nl 



6 {NcNcf ' 
NabN^^N^Nc . 

25 NabcN^N^N^ ^ 
(6/nl/5)2- 



(22) 

(23) 

(24) 
(25) 



B. General Single-Field Models 

For a single scalar field <f perturbing N{(p) we may 
expand using the above formalism [ISJ, as: 



(26) 



where N' = dN /dip. Note that we do not require that p is 
the inflaton field, it could be the curvaton or a field which 
modulates the efficiency of reheating. From equations [22]- 
24 we may immediately read off 



/nl 

3NL 



5 N" 



{N"f 

25 N'" 
54 (AT' 
1. 



\3 ' 



(27) 
(28) 

(29) 
(30) 



Biiki,k2,k3) = ^hL[Pdki)Pdk2) + 

Pc{k2)Pcik3) + Pdh)Pdki)], 



(20) 



T{ki,k2,k3,k4) = T^i^[Pdki3)Pdk3)Pdki) + in perms)] 
54 

+ ^5NL[Pc(fc2)Pc(^3)^c(fc4) + (3 perms)],(21) 

where Pdk) = NaNbC^^ {k) and therefore in the slow 
roll limit Pdk) = NAN^P{k). 

From the above two expressions we can read off the 



Equations [27] and [28] yield a very important conse- 
quence of single-field models namely tnl = (6/nl/5)^. 
This is a general result and therefore Anl 1 nray be 
used to rule out single-field models all together. 



C. Multi-Field Inflationary Models 

Suyama and Yamaguchi showed in general tnl ^ 
(6/nl/5)^ by the Cauchy-Schwartz inequality and equal- 
ity only holds if Na is an eigenmode of Nab 2Q] ■ Models 
where equality does not hold can not be those of a single- 
field. We now examine such models. 

Unlike the single-field case, using the SN formalism 
to make general statements about multi-field models is 
nearly impossible. Instead, one is forced to work with 
specific models that utilize simplifying assumptions. We 
now present a class of multi-field models that we believe 
is sufficiently general to uncover many details that are 
characteristic of multi-field models in general. 

Recently, Byrnes and Choi reviewed two field models 
with scalar fields ip and x that have a separable potential 
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W{ip,x) = U{ip)V{x) PinJH]. The slow roll parame- 
ters for these models are: 



the form 



~ 2 [ U 2 [v 



(31) 
(32) 



'Iw — ^'^p TT ' 'hx —"^p ' 'ixx —^"p y 



If ' -'^x —p ^ 

from which we can define 



V 



(33) 



For this class of models, in the regions where I/nlI > 1 
we have 



J, 5 f 

/NL = g^xTJ 



5NL = 

tnl = 

^NL = 



(l + f)2 

10 rjrivv - '^Vxx) - Vxx 
3 1 + f 

1 + r /6/nl' ' 



?~ V 5 

l+7~ 



r 



gHvy,.p-Vxx)^ ■ (34) 

/nl; (35) 
(36) 
(37) 



It is worth noting that both t^l and ^nl 8J:e related to 
/nl for this class of models. Here we have |(/nl| < |/nl| 
which will therefore be much harder to detect. On the 
contrary, tnl > (6/nl/5)^ so that non-Gaussinity may 
in fact be easier to detect in the trispectrum than the 
bispectrum for some multi-field models. Here we find 
^NL = {l + f)/f > 1. The scale dependance of /nl bas 
also been worked out for this class of models and was 
found to be Uf^^ = ~4(77^^ - r]^^)/{l + f) < 0. 



D. Curvaton Models 

In the curvaton scenario, a weakly interacting scalar 
field X exists in conjunction to the inflaton (p [6l [T8j [25l - 
[28] . During inflation, the curvaton field is subdominant, 
but after inflation x can dominate the energy density. 
The decay of the inhomogeneous curvation field in this 
scenario produces the curvature perturbations and not 
the inflaton. 

If such a curvaton field is the soul contributor to cur- 
vature perturbations, we can write out the perturbations 
using the 5N formalism as we did in the single field case: 

C = N'Sx + In"5x^ + iiV"'<5x' + (38) 

where now N' — dN/dx- Immediately we recover the 
relations [27]|29| and find for such curvaton models ^nl = 
1 as should be expected from curvature perturbations 
generated by a single-field. 

Recently, curvation models with generic potentials of 



(39) 



have been analyzed [27j [28] . Here m is the curvaton's 
mass and A is a coupling constant. For such models N 
in Equation |38] has been worked out giving: 



NL 



ffNL 



4r, 

25 
54 



-(l + /^) 



5?^ 
6 



4r^ 



(40) 
(41) 



+ -{l-9h) + 10r^ + 3rl 



where 



3n 



X,D 



XoXq 
X'o' ^ 



x'i 



(42) 



Here ^x,d is the energy density at time of curvaton decay, 
Xo is the curvation field during oscillations just before 
decay and the primes here denote derivatives with respect 
to time. 

Unlike single scalar field inflation, curvaton models can 
have large self interactions. Enqvist et al. pointed out 
that even if /nl is small, (7nl can be large for significant 
levels of self-interactions [2B: . This places a physical sig- 
nificance on gNL that can be thought of as parameterizing 
large self-interactions. 



E. Brief Summary 

In this section we have discussed the physical signif- 
icance of each statistic /nl, ffNL, tnl and Anl- In 
the bispectrum, /nl receives contributions from differ- 
ent shaped triangles in Fourier space related to different 
underlying physics. By analyzing the amount of non- 
Gaussianity from these different shapes we can distin- 
guish between models with multiple fields, non-canonical 
kinetic energy and non-adiabatic vacuums. 

In addition we stressed the physical significance of local 
non-Gaussianity in the trispectrum. The relation Anl = 
'''nl/(6/nl/5)^ is an important constraint of multi-field 
models. A general result for single-field models is Anl = 
1. Lastly, (7nl will place important constraints on the 
level of self-interactions. 



III. POWER SPECTRA ESTIMATORS FOR 
FIRST AND SECOND-ORDER 
NON-GAUSSIANITY. 

We would like to find a way to measure the non- 
Gaussianity of these fields from something directly ob- 
servable. Fortunately, information about the curvature 



6 



perturbations are contained within the CMB through 
the spherical harmonic coefficients of the temperature 
anisotropics: 



air. 



47r(-i)' 
ST 



(27r)3 



^k)gTi{k)Yr{k), (43) 



^(n) = — (n) = ^a,„yr(ft), 



(44) 



Im 



where $(k) are the primordial curvature perturbations, 
gTi is the radiation transfer function that gives the an- 
gular power spectrum C/ — {2/n) J k'^dkP^{k)g^i{k), 6 
is the field of temperature fluctuations in the CMB and 
are the spherical harmonics. (In this equation, the 
curvature perturbation $ is related to ( through the re- 
lation $ = (3/5)C.) 

If the curvature perturbations are purely Gaussian, 
all the statistical information we can say about them 
is contained in the two point correlation function 
(<I>(xi)$(x2)). The information contained in the two 
point function is usually extracted in spherical harmonic 
space, leading to the power spectrum C;, defined by: 



Ci 



1 



{21 



1 ^ 



O'lmO'iri 



(45) 



However, if the curvature perturbations are slightly 
non-Gaussian, this two point function is no longer suf- 
ficient to articulate all the information contained in the 
field. With non-Gaussianity, extra information can be 
extracted from the three, four and higher n-point corre- 
lation functions [5]- 

We now discuss estimators that can be used to measure 
non-Gaussianity at first and second order corresponding 
to the third and fourth-order in curvature perturbations 
respectively. 



A. Skewness Power Spectrum Estimator for the 
Bispectrum. 

In order to detect Gaussianity at first order, we must 
turn to the three point correlation function of the pri- 
mordial curvature perturbations ($(xi)<I>(x2)$(x3)). As 
mentioned above, we can extract information from the 
curvature perturbations by analyzing the ajmS of the 
CMB. The three point correlation function of the aimS> is 
called the bispectrum can be decomposed as follows [30]: 



{aimai'm'a,l"m") = B, 



ll'l" 



I V I" 

m m! ni" 



(46) 



where 



_ {2l + l){2V + l){2V + l) fll'l"\, 
Here the symbols in parenthesis are called the Wigner- 
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FIG. 3: The top plot compares various /3(r) for different r, 
and the bottom is the same for a{r). 



3j symbols and enforce rotational invariancc of the CMB, 
as well as ensuring the proper triangle equality holds be- 
tween I, V and I" namely: \li — lj \ < Ik ^ \h + lj\ for any 
combination of i, j and k. For more information on the 
wigner 3j symbols, the reader is directed to the appendix 
of Ref. [SOj. 

The quantity bwi", known as the reduced bispectrum, 
encases all the other information in the bispectrum and 
for the local model can be computed analytically as: 



2/. 



NL 



I 



where 



Mr) 



r^dr [ai^ ir)f3i.^ ir)(3i3 (r) + cyc.perm 



k'^dkgTi{k)ji{kr), 



k^dkP^{k)gTiik)ji{kr). (49) 



Here, P<s,{k) (x is the primordial power spectrum of 

curvature perturbations, griik) is defined above, ji(kr) 
are the spherical Bessel functions and r parameterizes 
the line of sight. 

Traditionally, the non-Gaussianity parameter /nl is 
given as a single number where the information from all 
the triangles configurations are collapsed into a single 
number called skewness (6*3) defined as 



j r'^dr j dhA{r, h)B'^{r, n). 



(50) 



where A{r,h) and B{r,n), are defined below in equa- 
tions [51] and [52] Recently, new techniques have been de- 
veloped to measure /nl from a power spectrum called the 
skewness power spectrum |3H I32j . These new estimators 
based on the analysis of power spectra are equivalent to 
53 in the limit of homogeneous noise [5] but have certain 
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advantages discussed at the bottom of this subsection. 
These advantages include the abihty to separate fore- 
grounds and other secondary non-Gaussian signals and 
the ability to correct for the cut-sky without having to 
compute so-called linear terms. 

To extract the skewness power spectrum from data we 
must begin with temperature maps optimally weighted 
for the detection of non-Gaussianity following [55] : 



Im 



Im 



Ci 

Mr) 
C, 



5;a;,„,(51) 



6;a/,„.(52) 



Here Ci = Cibf + Ni where bi and Ni are the beam trans- 
fer functions and noise power spectrum respectively as 
described below in Section |IVB| and C/ is the usual two 
point correlation function defined above in equation |45[ 

From the two above weighted maps we can create two 
unique two-one power spectra, each of which contribute 

(2 1) 

to the full C, ' estimator defined as: 



A.B'^ 



AB,B, 



(53) 
(54) 



where 



(r) 



C,^^'^(r) 



21 



21 



^5]Real{A,„,(r)i?,(^)(r)}; (55) 

m 

^ Real {iAB)imir)Bimir)}(,b6) 



It should make sense that the integrals with respect to 
the line of sight are needed since the final power spectra 
must only be an I dependent quantity. 

In the above equations, the sgMarerf multipole moments 
are defined in relation to the squared optimized temper- 
ature maps as: 



Bir,nf = Y.Bl'J{r)Yr{ny, 



(57) 



Ini 



Air,h)Bir,n) - ^(Ai?)„„(r)l^"*(n). (58) 

Im 

Combining the two unique contributions from equa- 



tions 53 and 54 gives us our full skewness power spectrum 
estimator: 



(2,1) 



A.B^ 



2a 



AB,B\ 



(59) 



(2 1) 

Once C; ' has been extracted from data, we can com- 
pute the amount of non-Gaussianity found therein by 
relating this estimator to its analytical expression for a 



model with /] 



NL 



1 that turns out to be: 



C, 



(2,1) 



NL 



{21 + 1) 



EE 

I' I" 



ii'i" 



Bu'i"B, 
CiCi'Ci" 



(60) 



Here, Ci is the weig hted two point power spectrum de- 
is the full bispectrum and 
1 calculated from 



52 



B, 



ii'i" 



fined below equation 
Biinn is the local model with /nl 
equations |47] and HSl 



Measuring non-Gaussianity using a power spectrum 
has a few advantages related to the fact that all infor- 
mation is not squeezed into a single number. First, dif- 
ferent physics that contribute to the bispectrum, such as 
point sources and secondaries, can be directly accounted 
for and measured using curve fitting techniques utilizing 
each quantities two-one spectrum and fitting all parame- 
ters simultaneously as was done recently in Smidt et. al 
2009 m. Second, each statistic can be tested for scale 
dependance with ease. This was also done in U where 
it was found that /nl is consistant with zero for all 
Third, effects due to the cut sky can be removed easily 
without needing to calculate linear terms needed with 6*3. 
We discuss this later issue in Section [IV C[ Lastly, for the 
trispectrum analysis discussed below, both second order 
statistics tnl and ^nl can be calculated simultaneously 
using the two kurtosis spectra. 



B. Kurtosis Power Spectrum Estimators for the 
Trispectrum 



In order to extract non-Gaussianity at second order we 
must consider the trispectrum or four point function of 
temperature anisotropics which conveniently breaks into 
a Gaussian and non-Gaussian or connected piece [34] : 

(61) 

(^/imi ^^2^2 '^^S^Tls ^^4m4 ) ^" ('^/l mi ^^2^2 ^^3^3 ^^*Tl4 ) ^ • 

where the connected and unconnected part of the trispec- 
trum can be expanded as: 

(62) 

y(-i)^^G;t(i) f t]( \r)^ 



(63) 

Ef-l)*^T'^'''fL) ( h L V I3, U L \ 
\ ) IJ2 ^ M mi 1712 M l\ ms 1714 —M I ' 
LM ^ ^ 
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where we can solve for G;'J;*(L) and Tj^^i^{L) analytically are: 



as: 



(64) 



/ 



rldrirldr2FL{ri , r2)ai^ ('"1)^2 {rijai^ {r2)l3ii {r2) 
r^drPi, (r) A, (r ) {r)Pi^ (r) +^^,3 (r) A, (r )] , 



jABA.B _ 


/ rjdri / rldr2Ji^^^'^{ri,r2); 


(70) 


^AB^,B _ 


f AB^ B 


(71) 


jAB.AB _ 


1 rjdri J rldr2jf'^''^^ {ri,r2); 


(72) 


^AB.B^ _ 


j r^drCf^'^ (r). 


(73) 



with tnl and ^nl being parameters of second order 
primordial non-Gaussianity (see discussion in next sub- 
section for more information). Written in this form, 
"^hh (^) a-bove is called the reduced trispectrum and con- 
tains all the physical information about non-Gaussian 
sources f5S|. The full trispectrum, in general, contains 
additional terms based on permutations of li. We approx- 
imate the full trispectrum with the reduced trispectrum 
since we will be optimizing the estimator with weights 
to measure a single term of the full trispectrum. There 
are additional cross terms in our analysis that we then 
ignore. The approximation we implement here is al- 
ready costly computationally and the lack of including 
extra cross terms associated with permutation, at most, 
causes our error bars on the non-Gaussian parameters 
to be overestimated. Furthermore, as we measure non- 
Gaussian parameters using the reduced trispectrum, we 
can direcly compare our results with the previous predic- 
tions that also utilized the same approximation |29[ I35j . 

In above, the quantity hi-^1,1, is defined such that 



(2/1 -1-1) (2/2 + 1) (2^3 + 1) f h h h 





47r 



and 



(66) 



Fi(ri,r2) = - / k^dkP^{k)jL{kr,)jL{kr2). (67) 



Here, Pg,{k) cx is the primordial power spectrum 

of curvature perturbations, the a{r), (3{r) and gTi{k) are 
defined above and ji{kr) are the spherical bessel func- 
tions and r parameterizes the line of sight. 

As with the bispectrum, we would like to figure out 
how to calculate power spectra that can be related to 
analytical expressions proportional to tnl and (7nl- To 
do this we begin with the same weighted maps defined in 
equations [51] and [52] which leads to the spectra: 

= (5/3)VNLJ,-'^^'^ + 2.gNL/:r''^, (68) 
Kl'''^ = (5/3) Vnl^.-^^'-'^ + 2gNL/:r (69) 
where the unique two-two and three-one power spectra 



Here Ji^^^'^ {ri,r2), Cf^ '^(r), J;'^^"^^(ri, r2), and 
£AB,B j-^-j g^^g angular power spectra of their respec- 
tive maps. For example Cf^ '^(^) is defined as: 

Cf '''''' (r) = ^ Y.(^B%^BL (74) 

m 

where {AB^)im andi?;*^ are defined analogously with 
equations [57] and [58] 

Once the kurtosis estimators have been extracted from 
temperature data, we can fit the two unknowns tnl and 
.9nl from the two estimators simultaneously by com- 
paring them to their analytical expressions with tnl = 
5nl = 1 that turn out to be |5S] : 



(2,2) 



{21 



-1 -1 rpl^l4 

+ l)4-(2/ + l) Ci, 



1 Tllll{l)Tllll{l) 



Ci,Ci,Ci, 



^(3,1) ^ 1 ^ 1 Tlll^{L)flll^{ L) 
' {21 + l)f-^ i2L + 1) CiAAsCi 



(75) 
.(76) 



where fl^l^H) is the full bispectrum and tI^[^{1) is the 
local model with tnl = (7nl = 1 calculated from equa- 
tion [Ml 



IV. FISHER BOUNDS 

A. The Ideal Experiment 

In order to determine the optimal error bars for these 
estimators we must properly calculate their signal-to- 
noise ratios. For the bispectrum, the signal-to-noise ratio 
takes on the simple form 



(2,1) 



J2{2l + l)C, 



{2,1) 



(77) 



(2 1). 

where Cj ' is defined above in eq 

For the trispectrum we must calculate the signal-to- 
noise for both /Cp'^' and IC\'^'^\ In a best case scenario, 
the two estimators above are not correlated. In this case 



60 
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the signal-to-noise for each estimator is: 



= Y.{21 + 1)/C 

(2,2) I 



(2,2). 



(3,1) ; 



(3,1) 



(78) 
(79) 



Given the positive definite nature of (S/N) , the signal- 
to-noise increases as one computes to higher I values. In 
fact, for the trispectrum it has been shown that {S/N}'^ ^ 
^max where Zmax represents the maximum I used in the 
analysis HHj- 

In addition to the estimators themselves being cor- 
related, contributions to the terms proportional to tnl 
and ^NL come from different quadratic contributions in 
Fourier space. This further allows us to calculate the 
signal-to-noise for each of these terms in each estimator 
by setting the other to zero. For example, we can deter- 
mine the optimal signal-to-noise for the tnl term from 

(2 2) • 

say the /Cj ' estimator by setting (Jnl — embedded in 
equation [?8j 

Once the signal-to-noise is known, we immediately 
have a bound on the optimal error bars for our estima- 
tors through the inverse square root. For example, if we 
wanted to know the optimal la error bar that can be 



placed on tnl from the /C; 
the Fisher bound as 

o'(rNL) = 



(2,2) 



estimator, we can compute 



(80) 



'(2,2) 



with the restriction on {S/N)^^^ 2) to tnl by setting ^nl 
in this calculation. 



B. The Realistic Experiment. 

In the above equations we assumed a "perfect" experi- 
ment with no noise or beam with a full sky. We now must 
take in to account that real world experiments have an 
inherent noise associated with the detector and a beam to 
characterize its angular resolution. Both the noise and 
the beam reduce the signal-to-noise. Furthermore, the 
mask yields a cut sky that must be dealt with properly. 

The noise is often reasonably approximated assuming a 
homogeneous spectrum calculating TV; from the following 
relation: 



"pix"pixi 



(81) 



where CTpix is the rms noise per pixel and ripix is the solid 
angle per pixel. 

For the noise calculation taking into the inhomoge- 
neous coverage of real world experiments and a cut sky 



WMAP 7-Year [V] 
Planck [143(GHz)] 
EPIC [150(GHz)] 




1000 1500 
I 



FIG. 4: Beam transfer functions. The frequency band used 
for each experiment is in brackets. 



Ni is to be calculated by: 



47r/sky A^obs(n) 



(82) 



where /sky is the fraction of sky observed and A^obs is the 
number of observations per pixel [3] . 

In addition to noise, realistic detectors have limits to 
their resolving power. The resolution limits of the instru- 
ment, encoded in the parameter 0fwhm which represents 
the full-width-half-max of the resolving power. We can 
map this information into harmonic space in the beam 
transfer function 6; 



exp (-^Vbc 
^'fwhm 



(83) 
(84) 



Beam transfer functions for the WMAP, Planck and 
EPIC experiments are plotted in Figure [4j As one would 
expect, a larger ^fwhm results in the suppression of in- 
formation on larger scales. 



Mission 


^FWHM 






Frequency 


Planck 
EPIC 


7.1' 

5.0' 


2.2 X 10-6 
8 X IQ-^ 


0.0349 
0.002 


143 (GHz) 
150 (GHz) 



TABLE I: Parameters used to calculate the simulated noise 
and beam transfer functions for the Planck and EPIC experi- 
ment 38, 39 . We obtained WMAP noise and beam function 
from publicly available data. 



Working in spherical harmonic space, it is easy to cor- 
rect our estimators /Cp'^-* and /C^'^' for the effects due to 
noise and the beam. All that must be done is to preform 
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FIG. 5: Noise and bi relation, Ni/bf, for each experiment plot- FIG. 6: The power spectrum Wi of the KQ75 mask, 

ted against C; taken from WMAP 7- Year best fit parameters. 
The frequency band used for each experiment is in brackets. 



the transformation: 



Cr 



(85) 



in the denominator of Eq. 75 and 75 As should be in- 



tuitively expected, a large amount of noise, or poor reso- 
lution will result in a smaller signal-to-noise. Therefore, 
how much the signal-to-noise is effected is related to the 
relationship between C; and Ni/bf. For Ni/bf » Ci, 
the signal is greatly diminished. The relation between Ci 
and Ni/bf for the WMAP, Planck and EPIC experiments 
is plotted in Figure [S] 



C. Mask And Cut Sky 

To remove cut sky effects using the traditional 5*3 esti- 
mator, many linear terms must be computed that must 
be subtracted off UHl ED- Furthermore, the num- 
ber of terms that must be computed grows for higher 
n-correlation functions. The difficulty arises because the 
cut sky effects are compressed into a single number, mak- 
ing it difficult to subtract out. 

One advantage of probing primordial non-Gaussianity 
with skewness power spectra is that we can use tech- 
niques pioneered by Hivon et al. to remove mask effects 
from the spectra [40j . This technique is relatively simple 
and works identically for correlation functions of arbi- 
trary order. 

When one uses realistic data, a mask M^(n) must be 
applied to an all sky map Mixi) to get rid of unwanted 
sources such as the galactic plane. This mask there- 
fore affects the aimS> derived from the all sky A{r, n) and 
B{r, fi) defined in equations 51 and 52 used in the bispec- 
trum and trispectrum analysis producing cut sky ci/mS: 



air. 



dhM{n)W{n)Yi"'*{h), 



(86) 



E 

I'm' 



ai' 



dhY!r{n)W{n)Yr*{n), (87) 



ai'm'Kiml'm'[W], 



I'm' 

Here ai'm' is for the full sky, M(n) represents an arbi- 
trary full sky map and Kimi'm' \W] now contains all the 
cut sky information. 

Hivon et al. showed that a power spectrum based on 
such masked data can be corrected by: 



(89) 



where Mu' is a matrix defined by 



97' I 1 . . 



I'' 



Here Wi is the power spectrum of the mask W^(n). The 
power spectrum for the KQ75 mask is plotted in Figure [6] 
and the corresponding Mu is plotted in Figure [7] 

Furthermore, it has been shown that any power spec- 
tra of rank c\^'''^ for any p and q can be corrected with 
the same method using Mui |36j . Thus, we can correct 
the skewness and kurtosis power spectrum estimators for 
the bispectrum (rank p = 2, q = 1) and the trispectrum 
(rank p ~ 2, q = 2 and rank p = i, q = 1) using this 
same technique. For example, a plot showing the effec- 
tiveness of this correction on the K^'^''^^ estimator is seen 
in figure [8] 

This correction technique is unique to the power spec- 
trum approach to detect non-Gaussianity because not all 
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Full Sky/Cut Sky Simulations with AIn, 






Cut Sky No Mif 




Full Sky 








BSSSisaaiiBiiiisg^ 



300 
I 

FIG. 8: On the bottom we see the results of 250 full sky 
gaussian simulations of IC\ ' with the cut sky results without 
correcting with Mui . On the top we see the same except the 
cut sky has received the proper correction. 



I dependent effects of the mask have collapsed into a sin- 
gle number. Therefore, the abihty to correct for the mask 
in this approach is much easier and more efficient than 
calculating linear terms needed to correct for masking 
effects in for the traditional skewness statistic 6*3. 



V. FISHER ANALYSIS AND RESULTS 

We now calculate the signal-to-noise for each of our es- 
timators in order to give reasonable expectations for non- 
Gaussianity detection from upcoming experiments using 
skewness and kurtosis power spectra. These constraints 
assume only temperature data from one frequency band 
per experiment is used. For the WMAP 7- Year analysis 



we use the V frequency band and for Planck and EPIC we 
use the 143 and 150 GHz frequency bands respectively. 
The noise and beam for the WMAP 7- Year V band was 
taken from the WMAP team and those for Planck and 
EPIC were computed using the values in Table |l] as de- 
scribed in Section IIVBI . 

It should be noted that combining different frequency 
bands and adding polarization can further reduce the ex- 
pected error. For example, with the recent WMAP 7- 
year findings error bars on /ml from one frequency band, 
V or W, is ±24 but the full temperature analysis com- 
bining V+W bands gives a reduced error bar of ±21. 
(About a 12.5% improvement over one temperature fre- 
quency band alone.) 

For each of these calculations Ci, a{r), /3{r) and 
FL{ri,r2) were calculated from Eq. 



49 and 



67 



usmg a 

modified version of CAMB based on the WM^P 7- Year 

best fit cosmological parameter values. The quantities 

Ci, a{r) and /3(r) are plotted in Figure [sj 

For the bispectrum we can form one skewness power 
(2 1) 

spectrum estimator ' which places bounds on the 
first order non-Gaussian parameter /nl- To calculate 
the signal-to- noise, we compute Eq. [77] from Eq.[60] sum- 
ming all I up to some Z^ax between 2 < Z < 1000. After 
calculating this signal-to-noise we calculate the Zmax de- 
pendent error bars from the Fisher matrix in eq. |80[ 

The results of this calculation are seen in Figjo] and 
shown in Table [III This calculation is done for the case 
of no noise nor beam, as well as with the noise and beam 
for the experiments WMAP 7- Year, Planck and EPIC. 
As, expected, the error bars drop for higher Z^ax until 
one reaches the limits of detection for each experiment. 
For the case with no noise, the error bars fall off as ~ 

For the trispectrum we can form two skewness power 



spectrum estimators, /Cp'^"* and /C 



-(3,1) 



For primordial 
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FIG. 10: On top we have Fisher bounds on tnl for the K\ ' 
estimator and on bottom for IC^ ' . The frequency band used 
for each experiment is in brackets. 
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FIG. 11: Top: Fisher bounds on gNL for the /Cf ' estimator. 
Bottom: Fisher bounds on (;nl for the IC] ' estimator. The 
frequency band used for each experiment is in brackets. 



non-Gaussianity detection that together place bounds 
on the second order non-Gaussian parameters tnl and 
(7nl. The first of these, /Cp'^\ is computed from eq. 75 
After this calculation, the signal-to-noise is computet 
from eq. 78 summing all I up to some /max between 
2 < I < 1000 for all I except the diagonal one. (The 
diagonal being the I in parenthesis of Ti\^{l).) It was 
confirmed, as was previously reported [29], that nearly 
all the signal-to-noise can be calculated only summing 
up the I in the diagonal of the trispectrum up to / = 10, 
saving a tremendous amount of computational time. In 
this analysis, however, we summed up the diagonal in 
both trispectrum estimators to / = 20 so as to be more 
conservative. The error bars on tnl and ^nl from this 
estimator are then computed from equation eq. |80l Re- 



sults from this estimator for tnl are seen in Fig. 10 As 

(2 1) . 

with the ' estimator above, we show the la bound 
for the case without noise and beam as well as for the 



WMAP 7- Year, Planck and EPIC experiments. For the 
case of no noise or beam, the error bars for this estimator 
fall off as - 1/V^nlP. 

Also plotted in the figure is the amplitude Anl assum- 
ing /nl = 32, the WMAP 7-year best fit value. There- 
fore, if /nl = 32 than we must have Anl > 1 for Planck 
to be able to have a detection of tnl- However, even 
if Anl ~ 1, EPIC should be able to detect tnl, espe- 
cially since EPIC will be able to use data much past 

/max = 1000. 

We also compute error bars for tnl from our second 
skewness power spectrum estimator for the trispectr um 

(3 1) . . . I I 

/Cj ' by first calculating the signal-to-noise from Eq. 176 



and 79 then solving for ct(tnl) from the fisher matrix [80^ 
Results for this calculation are plotted in Fig. [TO] Along 
with the la error bars for each experiment, is the am- 
plitude Anl assuming /nl = 50. The purpose of setting 
the amplitude to this value is to demonstrate that if /nl 
is large enough, models with Anl < 1 may be able to be 
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tested by upcoming experiments, especially EPIC. 

In addition to tni, bounds can be put on (/nl from the 
two before mentioned four point estimators. To do this, 
we calculate the estimators from eq. [75] and [76] setting 
tnl = and ^nl = 1- From here, we calculate the 
signal-to- noise from eq. |78| an d|79| whereupon we compute 
Fisher bounds fr'om equation |80[ The results are seen in 
Fig. [11] 



Combining the two estimators IC\ ' and iCf''^' gives 
the minimum error bars for tnl and ^nl seen in Table [IT] 
as well as Figure [T2j These are comparable to those 
of [5^ and [3S] who calculated Fisher bounds assum- 
ing only cosmic variance limited sky. They did not use 
the power skewness estimator, however, their estimator 
is equivalent in the limit of homogeneous noise [3] . Kogo 
and Komastu [21] found a higher signal-to-noise than did 
Okamoto and Hu [35]. This paper finds a signal-to- nose 
in between these values. 



^(34) 



^inax 


500 


1000 


1500 


2000 


/nl Planck 


16 


10 


8 


8 


EPIC 


15 


7.5 


5 


3 


Tnl Planck 


4350 


1640 


1550 


1550 


EPIC 


3700 


920 


400 


225 


^NL Planck 


1.6 X 10^ 


1.4 X lO'^ 


1.3 X 10^ 


1.3 X 10^ 


EPIC 


1.5 X 10^ 


1.1 X lO'^ 


8.4 X lO'' 


6.0 X 10^ 


^NL Planck 


3.0 


1.1 


1.0 


1.0 


EPIC 


2.5 


0.6 


0.3 


0.15 



TABLE II: The minimum error bars at la for /nl, tnl and 



gNL, using both K.f''^^ and K-Y'''' estimators, for the Planck 
and EPIC experiments up to /max = 2000. As stated in text, 
this assumes only one temperature frequency band is used in 
the analysis. 



(3.1) 



From this table we see that tnl can be detected at 
95% confidence level by Planck if tnl > 3000 and EPIC 
for Tnl > 600. If /nl — 32 in the bispectrum, this 
equivalently means tnl can be detected if Anl > 2 and 
^NL > 0.4 respectively, again alluding to the fact that 
EPIC will be able to test some inflationary models with 
^NL < 1- 

Furthermore, as can be seen in Fig. |13[ for large 
enough Anlj the trispectrum is more sensitive to non- 
Gaussianity, even for Planck. This may turn out to be 
very important as some models predict Anl > 1- It is 
therefore imperative that Planck examines the trispec- 
trum for non-Gaussianity as it may turn out to be more 
likely to get a detection there than in the bispectrum. 

Some models predict an undetectable amount of non- 
Gaussianity in the bispectrum (For example, /nl ^ 1) 
with a large amount of non-Gaussianity in the trispec- 
trum. These plots let us know just how hig Anl must be 
in order for a detection of non-Gaussianity to be made 
in the trispectrum for such scenarios. 

From these plots we see, for /nl = 1, the trispectrum 
becomes more sensitive to non-Gaussianity than the bis- 





FIG. 12: Fisher confidence intervals for /nl, 3nl and tnl. 
The dark and light blue represent the 68% and 95% intervals 
respectively for Planck. The red and orange represent the 
68% and 95% intervals respectively for EPIC. 



pectrum at I = 1450, 830, and 500 for ^nl = 50, 90 
and 120 respectively. For /nl = 32, the trispectrum has 
more sensitivity at I = 2350, 1150, and 450 for Anl = 1, 
3 and 10 respectively and for /nl = 50 we have more 
sensitivity at I = 1500, 750, and 300 for Anl = 1, 3 and 
10 respectively. 

1)/AAnl for both Planck and 



Figure 14 shows (Anl 
EPIC. In this plot it is clear that both Planck and EPIC 
are in a position to rule out single field infiation by deter- 
mining Anl 7^ 1- Large sections of the parameter space, 
consistent with current measurements, will rule out Anl 
equal to unity by > ba. 

Note from table [IT] that the expected bound on ^nl 
is about two orders of magnitude weaker than that on 
Tnl, even though both parameters are suppressed by a 
power spectrum cubed in (21 1. One reason is that the k 



dependent shape factor multiplying tnl in (21) diverges 
whenever kij — > 0, while the factor multiplying gNL only 
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FIG. 13: Comparison of the sensitivity of both the bispectrum 
and the trispectrum to non-Gaussianity assuming different 
values of /nl and Anl- 



diverges when one of the fc^ (and in this case the 
same applies for tnl as well). 



VI. PRIOR ANALYSIS USING THESE 
ESTIMATORS 

These skewness and kurtosis power spectrum esti- 
mators have recently been employed to constrain non- 
Gaussianity in the WMAP 5-year data. Using the bispec- 
trum, Smidt et al. (2009) found that -36.4 < /nl < 58.4 
at 95% confidence [4^ . This bound puts the la error bars 
at ±23.5, within about 12% of the optimal Fisher bound. 

The analysis for the trispectrum is more difficult and 
we therefore elaborate about it here. Our recipe for anal- 
ysis is 



(3,1) 



1. We calculate /C 
and ^NL = 1- 



2. We extract /cf'^^ 
5-year data. 



and /C 



(2,2) 



in Eq. 



75 



76 



for TNL 



and /Cf '^^ directly from WMAP 





FIG. 14: The top plot shows (tInl - 1)/AAnl for Planck and 
the bottom for EPIC. The color bands show to how many 
sigma j4nl would differ from unity for possible best fit values 
for Planck and EPIC. If Planck or EPIC find best fit /nl and 
Tnl values anywhere in the white region, single-field inflation 
will be ruled out by > 5cr. The black ellipse marks the 68% 
confidence region for the Smidt et al.(2009) best fit /nl and 
Smidt et al.(2010) best fit tnl values respectively g]. 



5. We fit the two unknowns tnl and ^nl from data 
using the two equations simultaneously. The am- 
plitudes the theoretical curves must be scaled by 
gives the values for tnl and (/nl 



3. We perform the extraction of /Cp'^"* and /C^'^'' from 
250 Gaussian maps, allowing us to determine error 
bars and the Gaussian piece of each estimator. 



6. We constrain j4nl by comparing tnl from the 
trispectrum with (6/nl/5)^ coming from the bis- 
pectrum. 



We subtract off the Gaussian contribution to these 
estimators to ensure we are fitting to the non- 
Gaussian contribution. 



This recipe is described in grater detail below: 

First we calculate /C^'^'' andjcj^'^^ theoretically using 
equations Eq. 64 ■ 67 and Eq. 75 - [76] for a model with 
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FIG. 15: The top plot shows the /Cf'^ and ICf'^ estimators, 
shown in green and blue respectively, taken from data for the 
W band. The same estimators for the V band are shown on 
the bottom. Additionally on the top the theoretical contri- 
butions for /Cf'^ and ACf'^ proportional to tnl are shown with 
the bottom showing those proportional to (?nl. The Gaussian 
contributions were not removed from these plots. 



TNL and ^NL = 1- To obtain d we use CAMB gT]^ 
with the WMAP 5-year best fit parameters and use the 
beam transfer functions from the WMAP team. We then 
obtain the connected piece using a modified version of the 
CMBFAST code gg^. Plots of many of the quantities 
used for these calculations can be found in Ref._p]. 



Plots of /cl^'^"* and JCf'^'' are shown in Figure 
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These 



curves will be compared with estimators derived from 
data to determine the magnitude of each statistic. Since 
we have two estimators, we can solve for the two un- 
knowns Tnl and ^nl by fitting both estimators simulta- 
neously. 

To calculate^ the estimators from data, used in the 
lefthand side of equations ( 75 1 and ( 76 1 , wc use both 



the raw and foreground-cleaned WMAP 5- Year Stokes I 
maps for V- and W-bands masked with the KQ75 mask ^. 
We use the Healpix library to analyze the maps. For this 
analysis we only considered data out to l^ax — 600. We 
correct for the KQ75 mask using a matrix Mui , based on 
the power spectrum of the mask, as described above. 



Figure |15| shows the results for /C^'^ and IC^'^ for the V 
and W frequency bands extracted from the raw WMAP 
5- Year maps. In order to do proper statistics for our data 
fitting we create 250 simulated Gaussian maps of each 
frequency band with riside = 512. To obtain Gaussian 
maps we run the synfast routine of Healpix with an in-file 



Gaussian Average 

Theoretical 

Data 



"iv.-- 



100 200 300 400 500 600 

I 



FIG. 16: The relation between the full estimators coming 
from data versus the Gaussian contributions. The green curve 
show the Gaussian contributions coming from averaging the 
estimators from the Gaussian maps. The red curve is the the- 
oretical Gaussian piece calculated using the WMAP-5 best-fit 
cosmology power spectrum. The error bars show two standard 
deviations from the Gaussian curves. These curves are from 
W band data. 



representing the WMAP 5-year best-fit CMB anisotropy 
power spectrum and generate maps with information out 
to I — 600. We then use anafast, without employing 
an iteration scheme, masking with the KQlb mask, to 
produce aim's for the Gaussian maps out to I — 600. 
Obtaining estimators from these Gaussian maps allows us 
to uncover the Gaussian contribution to each estimator in 
addition to providing us information needed to calculate 
the error bars on our results. 

This whole process is computationally intensive. To 
calculate all theoretical estimators took nearly 8,000 
CPU hours. Furthermore, all the estimators from Gaus- 
sian and data maps combined took an additional 1600 
CPU hours. 

As previously discussed, the full trispectrum can be de- 
composed into both a Gaussian and non-Gaussian or con- 
nected piece. To make a measurement of non-Gaussianity 
we to subtract off the Gaussian piece from the full trispec- 
trum. Figure [16] shows the the relationship between the 
full trispectrum and the Gaussian piece. In this plot the 
Gaussian piece was calculated in two different ways as 
a sanity check. First, the Gaussian maps were averaged 
over. Second, the Gaussian piece of each estimator is 



calculated theoretically using Eq. 65 



^ http : //camb.info/ 

^ http ; //www. cfa.harvard.edu/ mzaldarr/CMBFAST/cmbfast.html 
^ see Smidt el al. 2009 for a similar calculation using the bispec- 

trum for more details. [4] 
* http : / /lambda. gsfc.nasa.gov/ 



After obtaining the theory, data and simulated curves 
we use the best fitting procedure described in [3] where 
we minimize to fit tnl and ^nl simultaneously. Our 
results are listed in Table |III| We see that ijnl and tnl 
are consistent with zero with 95% confidence level ranges 
-7.4 < ^nl/IO^ < 8.2 and -0.6 < tnl/10* < 3.3 for 
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Band 


W 


V 


V±W 


Raw 








9NL 


4.7 X 10" ±5.3 X 10^ 


4.6x10* ± 5.9 X 10^ 


4.7x10* ± 3.9 X 10^ 


tnl 


(1.63 ± 1.27) X 10* 


(1.68 ± 1.31) X 10* 


(1.64 ±0.98) X 10* 


^NL 


7.4 ±7.3 


6.3 ±6.0 


11.1 ± 7.3 


FC 








Pnl 


4.2 X 10'' ±5.3 X 10^ 


4.1 X 10* ± 5.9 X 10^ 


4.2 X 10* ± 3.9 X lO-""' 


Tnl 


(1.32 ± 1.27) X 10* 


(1.39 ± 1.31) X 10* 


(1.35 ±0.98) X 10* 




6.0 ±6.7 


5.2 ±5.7 


9.2 ±6.1 



TABLE III: Results for each frequency band to la. Values for Qnl, tnl and tInl on the top are for raw maps. The values 
on the bottom are for foreground clean maps. Anl is estimated assuming /nl = 32 ± 21 from the WMAP-7 analysis and the 
tabulated la uncertainy is based on an analytical error propagation. 




FIG. 17: The 95% confidence levels for (/nl versus tnl. The 
red and orange represent the 68% and 95% intervals respec- 
tively for the combined V+W analysis. The light blue regions 
represent the 95% confidence intervals for the V band analy- 
sis, and the light green regions are for the W band. 



V+W-band in foreground-cleaned maps. The 95% confi- 



dence intervals of gjssL versus tnl are plotted in Figure 17 
for each band. For a V band analysis alone, the 68% 
confidence intervals are tnl = (1.39 ± 1.31) x 10'' and 
5nl = 4.6x 10'*±5.9x 10^. These error bars are within 
~40% and ^20% of the optimal Fisher values discussed 
above comparing with WMAP 7- year level noise for tnl 
and gNL respectively. 

Combining /nl — 11 ± 24 from Ref. [4 and tnl — 
(1.35 ± 0.98) X 10* from our skewness analysis we get 
—649 < Anl < 805 at 95% confidence. If instead we 
had assumed /nl = 32 ± 21 from WMAP-7 analysis [3] 
and same tnl reported here we find — 3 < Anl < 21.4 
at 95% confidence. The difference of the two estimates 
is a reflection on the central value of /nl since Anl = 
'''nl/(6/nl/5)^ and therefore a smaller /nl results in a 
larger uncertainty in Anl- This behavior is also seen in 
Fig. [Ml 

No measurements involving WMAP 7-year data have 



been preformed using these estimators. It is our opinion 
that the results for WMAP 7-year data will not be much 
different than for the WMAP 5-year data, just as the 
optimal results using the traditional skewness statistic S3 
do not differ significantly between these two data sets. [H 
i 

Planck, on the other hand, is in a position to make 
significant improvements in the measurement of non- 
Gaussianity using these estimators. Since Planck is tak- 
ing data, we encourage any plans to measure /nl, <?nl 
and Tnl using the skewness and kurtosis spectrum statis- 
tics that we have proposed. In addition to ruling out 
the standard single-field slow-roll inflation model with a 
detection of non-Gaussianity in general, Planck is in a 
position to possibly rule out all single-field models with 
a measurement of ^nl 7^ 1- 



VII. CONCLUSIONS 

In this paper we discussed the skewness and kur- 
tosis power spectrum approach to probing primordial 
non-Gaussianity. We outlined the expected constraints 
these techniques will place using future experimental 
data. These constraints were calculated by computing 
the signal-to-noise for each estimator, properly taking 
into account the noise and beam of each experiment. 
Optimal error bars for /nl, <7nl and tnl are listed as 
a function of /max- 

It was argued that the skewness and kurtosis power 
spectrum approach to measure non-Gaussianity has sev- 
eral advantages. These advantages include the ability to 
separate foregrounds and other secondary non-Gaussian 
signals, the ability to measure the scale dependance of 
each statistic and an advantage that the cut sky can be 
corrected from a matrix Mui without needing to compute 
extra linear terms. 

The physical significance of each non-Gaussian statistic 
is discussed. In the bispectrum, different non-Gaussian 
triangle configurations in Fourier space contributing to 
/nl are related to different underlying physics. By 
adding a local measurement of the trispectrum, a new 
statistic Anl = tnl/(6/nl/5)^ will be a powerful probe 
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to distinguish between multi- field models. Single- field 
models can be ruled out in general if Anl 7^ 1 and we 
discussed how this may be a real possibility with Planck 
or EPIC. Furthermore, for Anl large enough, the trispec- 
trum becomes a better probe for non-Gaussinity than 
the bispectrum for analysis utilizing information on very 
small scales. The parameter g'^L will be the hardest to 
constrain. A constraint on this parameter will uncover 
information on self-interactions. 
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